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Abstract 

We consider the spin n/2 analogue of the XXZ quantum spin chain. We compute the entanglement 
entropy S associated with splitting the infinite chain into two semi-infinite pieces. In the scaling 
limit, we find S ~ ^ ln(£) + ln(g) + • ■ • . Here £ is the correlation length and c K = is the central 

charge associated with the sfe WZW model at level k. ln(g) is the boundary entropy of the WZW 
model. Our result extends previous observations and suggests that this is a simple and perhaps 
£f~) , rather general way both of extracting the central charge of the ultraviolet CFT associated with the 

scaling limit of a solvable lattice model, and of matching lattice and CFT boundary conditions. 
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In this letter, we consider the spin n/2 analogue of the XXZ quantum spin chain (which corresponds 
to the choice k = 1). We derive an exact formula for the entanglement entropy S^ i,K ^ associated with 
splitting the infinite chain into two semi- infinite pieces. We consider the scaling limit and find that 
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s&k) ~ 2i\ n ^) +]n ^)) + c K . (i) 
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The first term involves the correlation length £ and c K = which is the central charge of the SI2 
WZW model at level k. The scaling limit of our quantum spin chain is known to correspond to a 
massive perturbation of this conformal field theory \ . This leading term involves the central charge 
in the way predicted in |2j and generalizes the k = 1 result derived there. The second term in Q is 
the 5I2 level k WZW model boundary entropy of Affleck and Ludwig 3 . It depends on the boundary 
conditions of our lattice model labelled by i € {0, 1, • • • , k}. The third term, C K , gives another finite, 
K-dependent contribution. We derive Q directly from the known expression for the partition function 
in terms of Ugfsl^) characters. 

The definition and meaning of entanglement entropy can be found in many places in the literature, 
see, for example, .4 . Our starting point is the density matrix defined by p = |vac)(vac| G End(7Y), 
where |vac) G H is the lowest energy eigenstate in the infinite tensor product space Ti of our quantum 
spin chain. The space TL is naturally a tensor product of semi- infinite left and right spaces Ti = Ji\®'H r . 
We are interested in the reduced density matrix defined as the partial trace 

Pl = Tv Hr (p) eEnd(Wi). 

The entanglement entropy we consider is then defined by 

S = -Tv ni ( Pl ]n Pl ). (2) 

In Baxter has discussed a range of 2D solvable lattice models which are related to ID quantum 
spin chains via a transfer matrix, and whose partition functions can be written in the form 

Z(x) = Tt Hi (x 2H c™), 
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where Hqtm, the corner-transfer-matrix Hamiltonian, acts on the left-hand space 7ii associated with 
the corresponding quantum spin chain, and has eigenvalues 0, 1, 2, • • • . It was shown in j^j, that for 
such models it is possible to identify the reduced density matrix pi discussed above as 



x 



Pi 



Z(x) 



A simple argument leading to this result is presented in .7]. The entanglement entropy defined by (J2J) 
is given by 
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S = -Tr Hl (piHpi)) = -^Tr Hl (pf)\ n=v 



which we can now write as 



S 



d 

dn 



Z(x n ) 
Z{x) n 



= ln(Z(x))-xln(x)|^ 

n=l Z\X) 



(3) 



This approach was used in ,2J in order to compute the leading contribution to the entanglement entropy 
for the Ising and XXZ models. In this paper we generalize these results to higher-spin models. 

The models of interest to us are the spin k/2 analogues of the 6- vertex model. These are solvable 
2D statistical-mechanical models defined in terms of spin variables e G {0, 1, • • • , k} that live on the 
edges of a square lattice. The Boltzmann weight of a spin configuration is given as the product of local 
weights R(Ci/C2) £ }' £ J associated with the following configuration around each vertex: 

ei 
Ci 



C2 



c~ 2 



The i?-matrix R(C) 1S understood in modern algebraic language to be an intertwiner of spin k/2 
evaluation representations of the algebra Uqfsh)- For explicit formulae for R(C) and a full discussion 
of this algebraic picture see [HI El - in the simplest k = 1 case, R(() gives the weights of the 6- vertex 
model, k = 2 corresponds to the 19-vertex model, etc. 

In this paper, we do not require explicit formulae for the Boltzmann weights, but just make use 
of known results for the partition function Z{x) (here x = —q, where q is the deformation parameter 
of the U q [5I2) vertex model). However, we are dealing with the infinite- volume model and do need 
to specify boundary conditions. We consider the model in the antiferromagnetic regime corresponding 
to the choice < x < 1, and consider a partition function, denoted Z^' K \x), which is defined to 
be the weighted sum over all spin configurations which are fixed at some finite but arbitrarily large 
distance from a nominal centre of our lattice to the following ground state configuration labelled by 
i G {0, 1, • ■ ■ ,k} (with i := k — i): 



The partition function Z^' L ' K \x) is given in terms of the principally specialised U q (sl2) character 
associated with the level k irreducible highest-weight representation V(Xi) with weight A« = iAi + (k — 
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z)Ao (where Aj and Ao are the fundamental U q (sl2) weights) |lUj . The explicit formula is 

Z(-)(X) = ^(A.)(^ = ^[f 2 ) +1)) - (4) 
Here, p = Ai + Ao is the f/^fs^) principal derivation (see 0), and we use the standard notation 

oo 

6 w (z) = {z]w)oo {wz~ l ] w)oo (w; w)^, and (z; w)^ = ^(1 - z w n ). (5) 

n=0 

A useful identity is 

(x , X )oo — (x , X )oo(^' , X ) QO (x , X )oo- (6) 

As an example, consider the k = 1 case. After using (jfifl. we obtain Z^x) = l/(x 2 ; x )oo- Note, that 
Z( 0,1 )(x) = Z^ l,l \x) in this case, reflecting the e = <-> e = 1 symmetry of the 6- vertex model. 
In this paper, we are interested in the entanglement entropy defined by 

^)=ln(zM(x))-xln(x)|^. (7) 

However, before computing this object, let us first discuss the scaling limit of the vertex models. The 
correlation length of the k = 1 model is given |12l |H] by 

where k' is the conjugate modulus associated with elliptic nome x (a gentle introduction to elliptic 
functions can be found in Chapter 15 of 5 ). If we define e by x = e _<E , then the £ — > oo scaling limit 
corresponds to e — > 0, x — > 1. In this limit fc' ~ 4e _7r ' e — ► 0, and so from (jHJ) we have £ _1 ~ A;' such 
that 

ln(C) k Ye (9) 

It follows from the results of ^3] that © also holds for general k. 

The series appearing if we substitute (@J into (JJJ) with x = e _e converges slowly for small e and is 
not appropriate for a consideration of the scaling limit. To proceed, we Poisson re-sum the series. Let 
us show how this works for a certain building block T(a, b) which appears in characters. We define 
T(a, b) for < a < b by 

T(a,b) := {x a -x h ) 00 {x b ' a -x b ) OQ . 
The contribution to the entanglement entropy from such a building block is 

S(a,b) = HT(a,b)) -xln(x) rfr( ; a ' 6) 

ax 

oo 

n=0 «ez 
where /(y) := ln(l - e v ) + _ . 

We can use the Poisson summation formula (see p. 468 of 5:) to rewrite S(a,b) as 

2 ^-/2vrn\ 2™a ~ f 00 vr / Try - sinh(vry) cosh(yn/)\ 
S M) = ^E/^rJ 6 6 ' where /(„):= ^ /(x) cos(yx)dx = - ^ j. 
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Noting that /(0) = 



— 7r 2 /3, we arrive at the expression 



n=l x 



The first term clearly diverges as e — > 0, whereas the series converges for e small. 

Using the fact that lim —f(S) = — vr/2, it follows that the limit of the series appearing in (I1UI) is 

5-^0+ 

given by 

* j b t ?(^) c- (T) - £ (^) - ^ (*- (")) • <"» 

n=l n=l 

Hence, in the scaling limit we have 

S(a,&)~-^ + ]*r2smRr)V (12) 



3eb V V 6 y 

with other contributions vanishing as e — > 0+. 

Now let us return to the computation of the entanglement entropy ([7)l. We can express 
entirely in terms of elementary blocks T(a, 6) by making use of © and the identification (x b / 2 ; x b )oo = 
T(b/2,b) 1 / 2 . We find 

ZM(X) = r r(2(i+l),2( W +2)) for ^ ^ 

T(2;4)2~ T(k+2,2(k+2))2 T(2,2(k+2)) T(4,2(k+2))-T(k,2(/c+2)) 

= 1 T(2(i+l),2( W +2)) for 

k odd. 

T(2;4)2 T(2,2(k+2))T(4,2(k+2))---T(k+1,2(k+2)) 

Each term T(a, 6) in the numerator will contribute S(a, b) to the entropy; each term T(a, 6) or T(a, 6)2 
in the denominator will contribute —S(a,b) or —-^S(a,b) respectively. Noting that all but one T(a,b) 
terms have the same value of b = 2(k + 2), and that a only appears in 1)10(1 via the cosine term, it should 
be clear that a simple exact expression for S^ 1 '^ then follows by using elementary combinatorics and 
trigonometric function identities. In the scaling limit, we use the form (|12|) and obtain 

S (i,n) = + ln V (*+ 2 ) / I for k both odd and even, (13) 

!2e \ (k + 2)2 

where c K = 3k/ (k + 2). 

The boundary contribution to the thermodynamic entropy of the SI2 level k WZW model was first 
obtained in [2] and is given by ln(</ l ' K )), where g( i>K \ the 'ground state degeneracy', is 

g (i, K ) _ sin (t^St 

Sill (l^2)) 

In the WZW context, the integer i labels the Cardy state of the boundary CFT [14] . It then follows 
that after making use of @ we arrive at expression JQ) for S^' K \ The term C K in Q includes the extra 
terms in the ln function in (|13j) as well as e independent, but possibly n dependent, corrections to ©• 
We would like to make several comments by way of conclusion: 

• The leading term in (^Q) is universal and involves the central charge in the way discussed in 
detail in [2^ I15| [TH] (the scaling limit of the spin-/c/2 model is known to correspond to a massive 
perturbation of the 5^ level n WZW model with central charge c K ). 
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• The second term in (^Q) depends upon the boundary conditions labelled by i. A similar contribution 
to the entanglement entropy from each boundary was found in [2] for critical model models on a 
finite length system with open boundaries (see equation (1.2) of j^j and try not to be confused 
by our reversion to the original ln(p) notation of |Hj - see also |17j). It is an interesting and 
apparently new observation that such a term is present for infinite size off-critical lattice models 
and survives in the scaling limit. 

• It is possible to derive Q by applying the conjugate modulus transformation x = e~ e — > x = 
e~ n / e directly to the theta functions appearing in equation (0J). We have chosen to present a 
rather more pedestrian approach in the hope that it may be applicable to a wider range of solvable 
lattice models, where the theta function realization of characters is not so well understood. 

• Our approach to computing the entanglement entropy should be generalizable to other familiar 
examples such as RSOS models, and hopefully to more exotic models. We hope that it might 
actually provide a simple alternative to the thermodynamic Bethe ansatz as a way of extracting 
the central charge of the ultraviolet CFT associated with a wide range of solvable lattice models. 
We also hope, given the presence of the boundary entropy term, that this technique for computing 
entanglement entropy might generally be useful for matching physical off-critical lattice boundary 
conditions with Cardy states in boundary CFT. 
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